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ABSTRACT 


This study describes two models for the computation of 
burst kill probabilities for air defense gun systems firing 
non-fragmenting projectiles at non-maneuvering aircraft tar- 
gets. Model I was suggested for U. S. Army use by Braddock, 
Dunn and McDonald, Inc. and is currently used in the TACOS II 
air defense battle simulation. Model II was developed by 
the Systems Analysis Directorate, HQ, U. S. Army Weapons Com- 
mand for use in gun system engineering development. The models 
are contrasted in development to demonstrate the strengths, 
weaknesses and relative merits of each. The Weapons Command 
model appeared to be based on a less restrictive set of 
assumptions than the BDM model, but sample results showed 
near equivalence in model BKP values throughout a reasonable 


range of engagement conditions. 
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I. INTRODUCTION TO THE PROBLEM 


Air defense gun systems have maintained an important role 
in air defense operations despite the advent of effective 
surface-to-air missiles. Reference 1 reports that the U. S. 
Air Force lost 550 aircraft to air defense guns in the Korean 
Conflict, and more recently, U. S. forces lost more than 3000 
fixed and rotary wing aircraft to enemy action in the Vietnam 
War with just a small percentage of these losses due to 
surface-to-air missiles. Furthermore, these losses do not 
reflect the considerable aircraft repair burden imposed by 
air defense guns nor the contributing effects (short of air- 
craft destruction) that guns have in deterring piloted air- 
craft from the accomplishment of their objectives. 

The air defense objectives have remained relatively con- 
stant over time, but the problems and technology of the gun/ 
target engagement have changed. Helicopters are now an inte- 
gral part of the battlefield, tactical aircraft are faster and 
more maneuverable, contour flying is possible under adverse 
conditions and air-to-surface weapon systems can be incredibly 
sophisticated. At the same time, digital computers, improved 
sensor capabilities and surface-to-air missiles which force 
higher targets to low altitude have favored the efficacy of 
the air defense gun system. 

Considering this continually improving capability of air- 
craft to deliver their ordinance effectively, it is imperative 


that every effort be made to assure that our forces are 
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equipped with the most effective air defense gun systems 
possible. 

The ultimate measure of success in the design and per- 
formance of an air defense gun system is the tested ability 
of the system to shoot down low flying hostile aircraft 
Within the combat environment. There exists, however, a 
more general need to have a measure of system effectiveness 
from a theoretical point of view that would allow for system 
evaluation and for comparative analysis of competing systems 
or sub-systems that could be applied in the development phase 
through operational fielding of the system. One such measure 
of effectiveness evolves from a system's ability to accurately 
deliver its ordnance to the target's point in space. This 
ability can be transformed into a kill probability based on 
target and gun characteristics which can be expressed as a 
Single shot or burst kill probability. 

It is this type of kill probability development as a 
function of gun system and target parameters that motivates 
this study. Computing the accuracy of a gun system and then 
confounding it with aircraft vulnerability data to generate 
a kill probability is a difficult theoretical and practical 
problem. There exists in the literature a profusion of for- 
mulas, tables and computational methods to assist in the solu- 
tion of special cases of this type of engagement, but there 
is a genuine need for unifying theory within which all 
special cases can be rooted. Also of interest is the increas- 


ing use of simulation to model combat and predict its results. 





Simulations require validated methods to compute air defense 
gun system kill probabilities based on prescribed battle sce- 
narios. 

The conclusion to be drawn from these comments is that 
the computation of system kill probability is a valuable tool 
which, as a current problem, merits immediate attention. It 
is the intent of this study to present and analyze two exist- 
ing models for the computation of gun burst kill probability 
(BKP) in order that the general and relative merits of each 
might be discerned. The two models were selected for comparison 
because they are both in use at the present time by different 
agencies within the Department of the Army. One is used as 
a sub-model in the air defense battle simulation model, 

TACOS II, and the other is used by U. S. Army Weapons Command 
agencies for analysis of air defense gun systems in engineer- 
ing development. It should be noted that this study represents 
the author's interpretation of these models and is not to be 
viewed as an official or final position report. 

Model I is a model suggested for U. S. Army use by Braddock, 
Dumm and McDonaid, Inc. 7 Rem 2/x It is a relatively basic 
model patterned as a two-dimensional model for salvo (or 
burst) firing of non-fragmenting projectiles. 

Model II is a model suggested for U. 5S. Army use by the 
Systems Analysis Directorate, HQ, U. S. Army Weapons Command, 
Rock Island, Illinois / Refs. 3 and 4_/. Its development is 
similar to that of Model I in that target vulnerability is 


represented in an identical manner in each and projectile 
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impact points are characterized as isolated random variables. 
The actual development of Model II seeks to improve on Model I 
by being based on a less restrictive set of assumptions. The 
separate and aggregate effects of these assumptions are pre- 
sented for comparison in this study. 

The general development of Model I is traced in Chapter II 
of this study with more detailed attention given to the signi- 
ficant assumptions required in the process of model formulation. 
This chapter identifies the theoretical basis for the general 
approach of the model, and comments are added to amplify those 
areas believed to be most critical to the model results. 

‘Chapter TII continues with the presentation of Model II 
and begins the comparative analysis of the two models. Con- 
trasts and similarities between models are presented, but 
with emphasis placed on the development of key aspects of 
Model II. Two eGnetnite cases are developed for this model 
after the general form of single shot kill probability (SSKP) 
has been shown. The two cases, A and B, contrast the effects 
of the assumptions regarding the magnitude of SSKP. Case A 
assumes that the product of the number of rounds in a burst 
(n) and the SSKP is small, while Case B assumes only the 
condition that SSKP is small. 

Chapter IV concludes the comparative analysis of the two 
models by stating possible advantages and disadvantages of 
each. The general impression created by the @ualitavre¢ 
analysis suggests that Model IIA is essentially the same as 


Model I, whereas Model IIB appears to be based on a less 


del 





restrictive set of assumptions than is required in either of 
the other developments. Quantitative results of one represent- 
ative test of the three model forms are presented in this 
chapter to illustrate the conditions under which Model I 

and Model IIA diverge. Somewhat surprisingly, the quantitative 
analysis indicates that although Model 1 and Model IIB are 
developed with two different approaches, the BKP results 

remain nearly identical over the broad range of parameter 


values exercised in the models. 
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II. BRADDOCK, DUNN AND MCDONAL INC. MODEL 
FOR BURST KILL PROBABILITY-MODEL I 


A. MODEL I INTRODUCTION 


The general approach in Model I is the development of 
Single shot kill probability for one gun system firing at 
one-maneuvering target which is then transformed into a burst 
kill probability in order that it apply more directly with 
the gun's normal mode of operation. This model was specifically 
developed to support preliminary analysis of air defense gun 
systems by examining the projectile/target relationships in 
the form of equations which describe the intercept environ- 
ment. Additionally, it is used as a sub-model in the TACOS II 
air defense battle simulation to describe gun BKP as a function 
of that intercept environment. 

The operation of the model involves the following steps. 
The gun system tracks its target, computes target character- 
istics, predicts a time-dependent intercept point in space, 
positions the gun and launches the projectile(s). The target 
is represented by a vulnerable area(A.) in a plane in space 
that is perpendicular to the slant range(R) between the gun 
and the center of the symmetrically represented vulnerable 
area. The center of the A, is the center of an (X,Y) coordi- 
nate system that represents the reference plane on which 
impact points for each projectile are measured and their dis- 
tributions described. Figure 1 illustrates the geometry Ot 


the gun/projectile/target relationships. 
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The remainder of this chapter describes the generation 
of the — by first developing SSKP and then transforming 
iat inte AF welFinc on the assumption of independence between 
rounds ies burst. The final form of the BKP equation is 
then reached through a discussion of the parameters which 


constitute the model. 


B. SINGLE SHOT KILL PROBABILITY 


The model assumes that projectile impact points in the 
target plane are normally distributed about each axis with 2 fue 
zero means and independence between the random variables Ko ( 
and Y which describe the projectile position in that plane. 

The assumption of normality relies on the Central Limit 
Theorem suggesting that error in the delivery of the projectile 
is the sum of a large number of independent error sources 

none of which contribute very much to the total error by 
themselves. Experimental data is required to validate this 
assumption. 

The fact that the respective distributional means are 
assumed to be zero has special implications when considering 
a model for measuring single shot effectiveness of a weapon 
firing a non-fragmenting projectile. A non-zero value of 
expected impact point could result in a case where higher hit 
probability results from larger projectile dispersion. Fig- 
ure 2 illustrates the point by depicting graphically the 


possible relationships among hit probability, mean impact 
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point and variance for a one-dimensional, normally distributed 
random variable. 

Cases 1 and 2 obviously result in greatest hit probability 
with smallest dispersion, but Case 3 demonstrates the possible 
existence of a Situation in which hit probability is maximized 
by making dispersion some finitely large value for a non-zero 
mean. This theory is easily expanded to two-dimensions with 
the same results. Reference 5 provides complete development 
of these cases as well as adding other insights into the effects 
on hit probabilities for a variety of mean, variance and 
distributional relationships. 

The normality assumption allows the representation of 


delivery error as a function of two probability density func- 


tions. 
1 (x = me 
fy(x) = eo. a 208 } (1) 
1 (y - 1 eo 
a —— = { - ——32 (2) 
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where Uy = Uy 


The independence of X and Y suggest the following joint 


density relationship: 


Py y(xsy) Py (x) P.(y) f fy (x) dx ° f fy(y) dy 


x oy 
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x a 
————— i  )  dkdy (3) 


where Ry and RY represent the range of target 


values for x and y respectively. 


This joint normal density representing projectile distri- 
bution (delivery error) in the target plane has utility in 
its present form, but it is transformed into polar coordinates 
since target vulnerability is represented in this model develop- 
ment by a circular area. It is also convenient to assume equal 
variances with respect to each axis in order that the joint OY | 
normal might be represented as a circular normal mre, 
Reference 6, however, suggests that the greatest variability 
in projectile error occurs along the azimuth(X)-axis where the |x °C 
target velocity component is greatest. This implies that the | 


variances are in fact not equatable. 


Continuing the transformation to polar coordinates: 
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The model assumes at this point that the target area seen ( 
from the gun may be approximated by a smaller vulnerable area Hl 
(A) which is represented geometrically as aed Tr cTe WakLiSuirs | 
total sumr™ewen) to the target vulnerable area. This assumption 


relieves a great mathematical burden in the development of 


BKP, but it creates another burden for the user who must 
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develop a reasonable technique for estimating A. Tor ua 
_——— a — ee 

given target type, aspect angle and range. The notion of 

vulnerable area is discussed further in Chapter III.B. 


Continuing the development: 


nega a | 
7 r such that r= A/a 


Equation (4) now becomes: 
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ee =o e.g | oePI 5} rarde (6) 
on O 0 0 0 
and finally, 
A 
SSKP = 1 - exp( - v-) (7) 
oT of 


C. BURST KILL PROBABILITY 


The extension of this model to account for BKP is facili- ; 
tated by the assumption of total independence among rounds lint | A 
4‘burst. It is not an intuitively appealing assumption. Guns ae 
with high rates of fire, as a minimum, suggest some correlation 
between rounds with the same aim point and separated in the 


firing sequence by only a very small increment of time. 
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Consider the simple case of a burst of two rounds. Let 
A and B represent the events that each round hits the target. 
The first round fired corresponds to A and the second round 


corresponds to B. General probability theory Suggests: 


P(AUB) = P(A) + P(B) - P(ANB) 
where P(A/]B) = P(A)-P(B/A) 


Under the assumption of independence between rounds: 


n 
Sa 
TRA py 


| 2 
h hy p deny 4} Recep p24 
P(B/A) = P(B) | sotalaay "keene 96 


- { p(pvB) > P(AUG) 4 


Le t rH, 
P(AUB) = P(A) + P(B) - P(A)+P(B) 7 
If the assumption is not valid, then the degree of corre- 
lation between A and B must be ascertained in order that its 
: ; Bee e, 
effects on the model be known. For example, if testing 4 (el) 2 
- ie a = } 


demonstrated that there was positive correlation between 

A and B, then the conditional probability of B given A would 
always be greater than the unconditional probability of B. 
The result would be a consistently higher estimate for BKP. 
The degree of this correlation would dictate the magnitude 

of the error in the general model caused by such an assumption. 
It remains possible, however, that the net effect of an 
erroneous assumption at this point might be negligible, so 
once again test data is required to support or refute it. 
Possible error notwithstanding, the development continues 
with n equal to the number of rounds per burst, each normally 


distributed and each having independent SSKP as developed. 
y td ra / I aT 
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P(all n miss) = (1 - SSKP)” (8) 


BKP = 1 - P(all n miss) = 1 - (1-SSKP)" 


~ [{ - a ee 4 
ies | lS el eee” 
from Eq. (7), bet IU a hw iF 
<)> 9 2e 
" n AL 
BKP = 1 - exp ( - ——~— ) (9) 
27 Se / 


This model asserts that BKP increases with an increase 
in target AL and with an increase in the number of rounds in 
a burst. The variance term suggests BKP decreases if the gun 
system is less accurate in its delivery capability. In its 
current form, BKP is a function of three parameters none of 
which causes counterintuitive impact on the value of BKP. One 
difficulty, however, arises when the user tries to incorporate 
this model into his scenario. Where does one acquire the 
appropriate values for n, A. and geo 

Burst size can be fixed at a reasonable level based on 
known firing rates or it may be simply varied in the equation 
for BKP and in field tests to determine the degree of its 
effect on BKP. Vulnerable area and variance are not so easily 
deduced. AY is not necessarily a function simply of target 
total area, but rather it is more probably a function of an 
aggregation of many factors relating directly to the target, 
and perhaps others attributed to the gun and its projectile 
characteristics. This issue is discussed further in the 


development of Model II in Chapter II1.B. 
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The variance in a gun system's delivery capability is 
another factor that has considerable impact on model results, 
but is difficult to determine quantitatively. Expressions for 
variance have been advanced with varying degrees of resolution, 
and they all share the commonality of being just estimates. 
This model considers only the aggregated dispersion due to 
aiming and ballistic errors rather than an in depth analysis 
of the many system functions that contribute to the net effect. 
The model form allows for easy adjustment as techniques improve 
for developing system variation components. Variance as 
represented in Eq. (9) is a function of slant range (R) bet- 
ween gun and target as well as being related to aiming and 
ballistic errors. Model form requires that variance be expressed 
in units such as square meters. If the aiming and ballistic 
contributions to variance are developed in angular units, 
aad if it is assumed that angular dispersion is constant over 
range, then Eq. (9) is easily modified to include the range 
Se Zs L$ o~ 


L 
ad 
me 


a 
- AA eer at 


factor. : eo ver 


> “(meters*) = Re (meters ’) ° g é (radians) 


such that 3 
= 
= (10) 





BKP = 1 - exp ( - x3 ) 


Consider the sources of variance in the projectile distri- 
bution. This model suggests, as does Helgert in Ref. 7, that 


variance may be thought of as an aggregation of independent 


variances due to ballistic dispersion and due to aiming errors. 
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Let oe and on represent these variance terms in square 
meters. The aiming error variance represents all those fac- 
ters Which contribute to the condition that the gun is not 
actually pointing exactly at the desired aim point at the time 
of firing. It constitutes effects due to components that com- 
pute target parameters, system servo mechanisms, alignment 
problems, gunner effects and others. These factors, as 

was discussed in a previous section of this study, are assumed 


to result in a normally distributed aiming error with zero 


mean and the variance term as just described. In general, 


oF (meters*) a R¢ (meters*) Cee (radians) 


Ballistic dispersion is a function of many projectile 
characteristics such as muzzle velocity and weight plus effects 
due to ambient conditions. Experience has supported the 
assumption that ballistic error be represented as a random 
variable with a normal distribution of zero mean and variance 


as given. In general, 


os (meters<) = R“(meters“) Oo sf (paeee) 


o 


Considering that the sum of two independent normal random 
variables is again a normally distributed random variable with 


parameters also additive, it is possible to conclude that: 


; eZ ; Z 
Oo ; (radians“) =e) o r og BR (radians ) 
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When applied to Eq. (10), the result is the final form 


fer Model I. i, 
; : Ue 
Me AL KN pyri > 
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III. U.S. ARMY WEAPONS COMMAND MODEL FOR BURST 
KILL PROBABILITY-MODEL II 


A. MODEL II INTRODUCTION 


Model II's development parallels that of Model I in that 
it first develops an expression for single shot kill probability 
and then transforms that into an appropriate form for burst 
kill probability. It constitutes an expansion of the Analytic 
Gun Model developed by the University of Michigan, System 
Research Laboratory for the U. S. Army Weapons Commnd. See 
Ref. 8 for the original development. 

The model asserts that target vulnerability can be repre- 
sented in a plane perpendicular to the slant range from the 
gun to the target at the point of predicted intercept. There 
is existing abadoiloey that would allow the representation 
of three-dimensional targets, but its adaptation at this point 
in the model would be unnecessarily cumbersome to its develop- 
ment. Reference 8 describes the methodology as well as amplifies 
the notion of vulnerable area. 

The projectile distribution in this target plane has a 
probability density function defined as fy y (x,y) where the 
x and y values relate to coordinate axes centered at the target 
center of vulnerability. Some analyses refer touthismommsin 
as intended aim point. This density describes the position of 


the round at the predicted time of intercept(t). 
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Fart 8 of this»chapter initiates the theoretical basis 
for this model by describing the development of single shot 
kill probability. The distributional form of the random 
variables that describe rounds impacting the target plane 
1s proposed and the notion of vulnerable area is addressed 
in detail. Once SSKP is presented, the concept of a burst 
center random variable is introduced. This concept is one 
method for improving upon the undesirable assumption of 
independence between rounds in a burst used in Model I to 
transform SSKP to BKP. It is at this stage of the develop- 
ment of BKP in Part C of this chapter that Model II is 
divided into two cases. The two cases, A and B, are necessary 
in order that a key assumption regarding the magnitude of 
SSKP might be approached from two different aspects. The 
development of each case is related with some final comments 
concerning seven neve in the theoretical bases for Model I 


and the two cases of Model_Il. 
B. SINGLE SHOT KILL PROBABILITY 


The general development of SSKP relies on several key 


assumptions. The first of which is that target vulnerable 


oe 


area remains constant over the time of the engagement, i.e. 


— 
ites 


A(t) = AL. Tenestyemn seaeery docs not allow this to be 
technically true, but the validity of BKP as a good measure 

of system effectiveness is not necessarily negated by its 
acceptance. Assume further that the total aircraft Silhouetted 


area (A,) can be systematically transformed into a representative 
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vulnerable area, and that one hit in that area produces a 
target kill. If the distribution of the projectile in the 
target plane has a large dispersion with respect to that 


area, then the small target approximation may be invoked 


me conclude tnat: 


SSKP = fy yxy) : A. 


where x and y are evaluated at target center (0,0). 


The small target approximation is noted here even though 
the necessary distributional assumptions for the projectiles 
have not been advanced. It is known that this model will 
assume a bivariate normal distribution for the projectiles 
in the target plane with zero means and no correlation. The 
result is a distribution of the form hypothesized in Model I 
and represented by Eq. (3). If the variances are large 
relative to the x and y values for target dimensions, then 
the quotients in the power of the exponential are approximated 
by zero and the double integral reduces to: 

ee ; 


“xy me 
< t®o 


ao 

where x+y = target area for a rectangular target which 
equates to A. The quotient in this approximation corresponds 
to the density function in the SSKP formula when it can be 


assumed normal with zero means, no correlation and evaluated 


at thesaime poant, (OmO). 
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Considering this approximation as an analog to an argument 
ee conditional probability may help to clarify its basis 
While at the same time providing some insight to the notion 


of vulnerable area. In general, 


P@i11) =—"P@ait) © P(kill7/nit) 


which when applied to this example suggests: 


A 
P(kill) = P(hit on A,) + —~ 


Ae 





where the quotient here is essentially a lethality function 
in which the assumption that AL is uniformly distributed over 
A. “S implicit. The quantification of AL is a potential 
problem for the user of the model. Ballistic Research Lab- 
oratories are said to provide data of this nature, but in 
order to effectively exercise the model, one must be aware 
of the subtleties involved in defining A. The analogy to 
conditional probability suggests that: 

P(hit on A.) = fy yxy) A, 


where it can be seen that the density function is the probabil- 
ity of hitting a vulnerable area of unit size. This factor 
is based exclusively on gun system characteristics and leads 
to a slightly modified perspective of the original equation 


for SenP. Now, 





A 
Vv 
SSKP = fy y(xsy) 2 A. K,, 
somthade Gn thesoriginal form for SSsKR, AL equates to the 


product of Aa an@ the ratio,of AL to A... This term provides 


Zo 





input for the effects due to target parameters as well as 
some consideration for projectile characteristics. 

The conclusion to be drawn from this is that SSKP in Bist A 
its present form may be computed by describing during the -- a. Az 


Ax 
course of an engagement the presented target area, the 
jm OTA, 


lethality relationship for the given gun and target and 
cae distribution of ails 3 ail the target A. . 

The user is still somewhat burdened at this point by 
the need to determine a reasonable lethality relationship. 
Model I assumes the simplest case of uniformly distributed 
vulnerability over total target area. If the target is 
depicted geometrically as a shoebox in space, then the gun 
system views a silhouetted area made up of three side-related 
components of the box. The target vulnerable area is then 
just a scaled down representation of that total area as 
illustrated in Figure 3. The mathematics of the develop- 
ment is facilitated by the conversion of A. to a circular 
area of equivalent magnitude. Parry, in Ref. 9, suggests 
a higher resolution methodology for determining AL as a 
function of gun/target aspect angle and total presented 
area. It assumes pockets with varying degrees of vulner- 
ability within the target profile as illustrated in Figure 4, 
The effect of a round impacting in the cockpit or engine 
compartment might produce more catastrophic consequences 
than a round impacting in an avionics section or some other 


less critical area. The notion of compartments of vulner- 


ability that may be aggregated to a single vulnerability 
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factor has its merits in theory and in practice may be 


applied to this model just as simply as less detailed tech- 
niques. 

Having considered the notion of vulnerable area, it 
remains for this chapter to develop an acceptable form for 
the density describing the projectile distribution in the 
target plane before the final form for SSKP can be presented. 

Consider the locations in three-space (R5) of both the 
target and the projectile at the predicted time of intercept. 
A three-dimensional coordinate system is proposed in Ry with 
its origin at the target center of vulnerability with axes 


™ YY, and 2. 


X = horizontal (azimuth) coordinate in the plane of the 
target 
a —>—] 
Y = vertical (elevation) coordinate in the plane of the 
target 
Z = range coordinate perpendicular to the plane of the 
target along the slant range from the target to the 
gun. 
It is assumed at the time of predicted amecmae PER v) gets 
(i) There exists an error vector (AA, AE,A R)(t) 
corresponding to the distance between the target 
center of vulnerability and the projectile. 
(ii) Each component of the error vector is the sum 
of random errors due to independent random 


n 
variables, e.g. 4A = 2% ( AA)i. The Central 
i=l 


BZ 





Limit Theorem is invoked such that the error 
vector may be assumed to have a multivariate 
normal distribution of dimension three and zero 


covariances. 


This is a major departure from conditions described in Model 
Sim that due towange error, 4R, the projectile will pro- 
bably not be coincident with the target plane at (t). Conse- 
quently, the error distribution in Rs at (t) must be projected 
into the actual target plane Ro. The projectile distribution 
an. Re corresponds to the multivariate normal distribution in 
Ry at a time when range-to-target and range-to-projectile are 
equivalent for the gun. Define a new set of coordinate axes 
in Ro as (X°',Y') which correspond to the X and Y axes in Re» 
respectively. 

At the time of predicted intercept, the target is described 
by the velocity components (i, oR, E.R, R)(t), and the projectile 


velocity components are (0, 0, Veet)”, where: 


> 
II 


target velocity in the X-direction in radians/second 
E = target velocity in the Y-direction in radians/second 
R = target velocity in the Z-direction in meters/second 
R = slant range from gun to target in meters 
Mss = projectile velocity in meters/second 


The non-maneuvering target is assumed to have a constant 
velocity vector over the time of flight, tos of the projectiles. 
The degree of effect of this assumption can be substantial 


if the target is changing velocity during the engagement, 


© 





but there is no effect on the model as long as it is under- 
stood that this model is restricted to the non-maneuvering, 
constant velocity target case. Theory has been advanced for 
the maneuvering target case in Ref. 4, but it will not be 
addressed in this study. 

The time between conditions of equal range and time of 
predicted intercept is represented by the random variable, aq t. 
It 1s easily assumed that At is small such that the projectile 


velocity, V_, may be considered constant over at, and the 


Dp 
relative velocity of the projectile with respect to the 
target, (R(t) + V(t), may be approximated by V(t) - = : 
Therefore, 


Returning to the error vector ( M, AE,AR)(t), the multi- 


variate normal assumption implies that: 


2 

A A ~ N(M,, oO = 
A E % N(M » O = 
o J 

A R N(M oO ) 
Z Z 


These random variables are projected into R, at t + ae 
This projection results in the creation of two new random 
variables 4A’ and A E' which relate azimuth and elevation 
errors to the X' and Y' axes, respectively. Considering 


Eq.(13), the error terms become: 
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iH 


AA (t +A) (Wes = \ to. A Rec) 


(jh + AB + a + R)(t) (14) 


co 


(AE -at’® E- RIC +) 


Am (t+ #4 €) 


(jE + a “A + R)(+) (15) 


jA' and ,E' are the sums of the independent, normal random 
variables~A, aE, and aR. Therefore, they have a bivariate 
normal distribution with possible correlation and with para- 


meters given by: 


aan N(ML", 9 Ss) 
a 
AE a N(M_ oy') 


The result of combining these two distributions is that at 


=A t, 
ili i 
UC) el dll oe =O cre co 
xX »Y aT Txt Tyr =) 4 Gly sg) 
(x-m_,1) * (x-M,0) (y-M +) 
{ age aa aire G+" 
mM y 
(y-M_+)* 
+ —— } } (16) 
Oo t 
J 


This function must be evaluated at some point in the 
target plane in order that it represent the engagement at 


intercept. It is assumed that the target is ‘perfectly' 


35 





located and that the center of vulnerability is in fact the 
aim point so that the projectile distribution is projected 
into the target plane with the effect measured over target 
center. Therefore, the function is evaluated at (x,y) = (0,0). 
When combined with the vulnerable area term, this yields the 


final form for SSKP; 





eS 
Sacer = AL ° fyty '(0,0) 
A. i M, @ 
a Oe aD ) 

27 O10 y' ¥1-02 -2(1-p*) x 
Nee M g M 8 ec 

- 20 + = ) (17) 
Sage, 7S y O y } } 


C. BURST KILL PROBABILITY 


Model I assumed at this point in its development that 
reunds in a ne are independent of each other in effects 
on the target. As was discussed, this is counterintuitive, 
but it has merit in that it conveniently eliminated some 
Pee dated obedaeles in the model development. Model II 
suggests that the negative aspects of this independence 
assumption may be lessened by considering an alternative 
solution to the transformation from SSKP to BKP. 

The alternative proposes that groupings of rounds that 
are either distinctly separated bursts or arbitrarily divided 


bursts be considered for collective effect. The centers of 


the bursts are defined as (R,,R5) and are assumed to be 


Se 
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independent, normally distributed random variables. The 
advantage to be gained from the introduction of this random 
variable is that the Central Limit Theorem confirms that the 
distribution of the burst-centers is_in fact normal regardless 
of the underlying distribution of the individual projectiles. 
By assuming the mean value of the projectile distribution 
to be the burst center, BKP may be developed by conditioning | 
on the known burst center distribution. 
The following terms are defined to assist in the develop- 


ment of the BKP equation: 


SSMP = Single Shot Miss Probability = 1 - SSKP 
BMP = Burst Miss Probability = 1 - BKP 
n = number of rounds per burst 
(R,,R,) = burst center coordinates in the (X', Y') plane 


(X,,X.) = position of the individual projectiles impacting 


inewi® ,Y‘) planet 


Assume that: 


(i) ee and X,, are independently distributed in the 
(FY plane. This facilitates bypassing the 
problems associated with trying to quantify the 
correlation between the two random variables. The 
mathematics which follow are also more tractable. 


Once again, it is difficult to assess the ultimate 





Innis ig identical with the (X,Y) notation used earlier. 
It is changed here to facilitate the matrix notation that 
will follow in the solution. 
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effect on the model results, but it is clear that 
increasing the number of assumptions of this type is 
narrowing the margin of difference between the develop- 
ments of these two models. This particular assumption 


implies: 


2 i 
Xs vN(Rs 40 i ) i= 22. ory dO) rd, 


(ii) Burst centers (R,,R,) are normally distributed with: 


2 
Ry vN(M,,0 3 ) Livsh 


R, vN(Myso & ) 


with correlation coefficient, p». 


2 
a anda > are dispérsion factors in the distribution 


of individual rounds due to ballistic errors and gun dynamics 
errors. The g : and g i variance terms are a measure of 
the dispersion of the burst centers in the target plane. 

The subscripts 1 and 2 replace x’ and y’ for convenience. 
Consider the development of BKP conditioned on the distribu- 


tion of the burst centers. 


BMP = E(BMP/(R,,R,)) (18) 


E( (SSMP)"/(R,+R,) ) 


il 


E( (1-SSKP)"/(R,,R5) ) (19) 
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| 
it neSskr is small, then (1-SSKP)” may be approximated 


ee ee 





Se . ene EE ee RL 


by the first two terms of its binomial expansion. This assump- 
tion essentially says that the gun system being considered 


is functionally incapable of performing its intended mission. 


ad 


For gun systems with high rates ee experience suggests 
that SSKP may be assumed small, but not n-SSKP. Development 
of the model under both of these conditional assumptions are 
traced here to illustrate the impact of each. 


First, if n-SSKP is small, Eq. (19) may be expressed as: 





BMP = E {(1-n-SSKP)/(R,,R,)} 
A RS R« 
ne 
= a {i- Z ’ exp {- — - 25) [in .R) } 
21 05 C5 205 Zz J. . 
Therefore, 
2 
mA Ry RS 
BKP = E{ -—--—-—— exp { - —s - 5)/ (RR, )) (20) 
27 0, 55 2 07 z 05 


The solution of this expectation is facilitated by the use 


et matrix notation. Let: 


2 

a ee 3 On 
V = 

2 

. ok oh 

R. | M 
R=|1 and M=| 

ee My, 


Consider the following solution to a general conditional 
expectation as a model for the form of BKP as it is trans- 


formed from Eq. (20) to Eq. (21). From Ref. 11, 
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E{X}=_f E(X/v=y} + fy(y) ay 


Therefore, 
oo nA 0 
BK P sae —_ - exp (- > R 5 il R ): 
on oT 0. oO 1 ae 
0 

1 a ) 1 

: ’ = = V-~(R-M); aR, dR 
ae ee exp ( een R-M); 

2n| v |* é iL Z 


ez) 
Banash, in Ref. 3, develops the solution to this double 


integral by manipulating its terms until it can be observed 
that it is the integral form of the moment generating function 


om a multivariate normal distribution such that: 


n-A I ‘aan i 
BKP = ——“——__ * exp (- v A ( J[-(I+A) tayy') 
2G Gis $ = 
2% val ee 
where I = 1 0 
OD 
Ze 
Z 3 1 
2 
O o/ %% 


and Vv. =0, if it is assumed that burst center means are zero. 


ewe result is that: 








n‘A 
BKP = v ‘ 
a2 [+A 2 (23) 
2 
a + 2 
where |I+A 2 Sie °3 /o, 0 Vie 
Zz 
0 cr Oy Os 
2 2 Zo 2 $ 
rc a). - (Chee: 
Oy. oD p 2 
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and the final form for Model IIA, 
re e 


nes, U 
| 2n( Castors (oy +05) )3 is 
\ 


Now consider the second case where only SSKP_is-considered — 





small. The development proceeds essentially the same as that 
of the first case, by conditioning on the location of burst 


centers. 
BMP = E( (1-SSKP)"/(R,,R,) ) 
Consider the expansion of (1-SSKP)". 


(1-SSKP)™ = 1-n:SSKP + +5 -SSKP* 


_ n(n-1)(n=2) nop pine? ‘SSKP? + ---- 


n 
= 2 


(3) (-SSKP)* 
J 


0 


Then, 


ri ‘ 
BMP = E( 2 (5) (-ssKP)4/ (R,,R,)} 


EB ogg (3) (-SSKP)I/(R, RQ) } 


j=0 


j 
~_—— 


| 
dlr 


= n LD (-1)9 { ItA. 
so (3) 








2m 07 O95 


exp (-4 y'A,(I-(I+A,)7A;) v ) 


41 





Banash shows in Ref. 4 that B = 7°, suee that: 





BKP = 1-BMP 
A J 
_ Ilo ee Vv 2 a | ees 
at “hha ca ssva| 
-axp(-2-§ y'A(I-j(I+4A) A) v_) (25) 


BKP is difficult to quantify in this form, however, it may 
be greatly simplified if the burst means once again are 
assumed zero and if the correlation between Ry and Ro is 
small enough to assume to be zero. 


The result of combining these assumptions with the general 


BKP formula is: 














J 
n . A A 
ee Be p(n) I) { ee ee 2 ( 26) 
( J 
in 97 O_O 
2 
where 
-$ 2a 
I+jA | 72 = —~+ £4... 1 
((o2+ joD(oF+ jor yp 


Adding this result to Eq. (26) produces the final form of 





Model IIB: | 
f Natt Ay >/7! ¢) A j 2 s 
| +1 V i 2 
| 270595 {Co 7 +5 0 3)o tj oy 


(27) 
When just one round is fired, this result for BKP is 
exactly the same as the result in Eq. (24) when it was 


assumed that n-SSKP was small rather than just SSKP. This 
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adds some degree of confidence to the mathematical manipula- 
tions involved in reaching the two final forms of the BKP 
equations. 

It is interesting to note at this point that although 
the developments of Models I and II have appeared to follow 
completely different tacks, it seems that Model I and Model IIA 
have produced the same result. It is perhaps more clear if 


one considers the following. In general, for small a, 


ee 


= 1-(1-2) = 2 
Therefore, from Model I, considering the original assumption 
that dispersion was much larger than vulnerable area, Eq. (11) 
may be restated as: 
nA 
ee ae er (28) 
a fe coe!) 

This result may be viewed as nearly identical with Model IIA 
as expressed in Eq. (24). The only possible difference being 
the manner of generating and expressing the variance components 
of projectile delivery. Thus, it might be deduced that Model I 
is implicitly assuming in its development that nx SSKP is 
small, or conversely, Model IIA may have implicit within its 
development the assumption of independence between rounds of 
a burst. The negative significance of these two basic assump- 
tions cast serious doubt on the veracity of each of these 


models as they exist in Eq. (11) and Eq. (24). 


a 





Nonetheless, the extension of Model II resulting in 
Eq. (27) provides alternatives to these two undesirable 
assumptions. Stating optimism for this modification to the 
model is one thing, validating its results remains another. 
Some discussion has been presented regarding the notion of 
vulnerable area since it is a key input parameter to the 
models, but very little information has surfaced in this 
model about the origins of the variance terms. The Appendix 
is included to demonstrate the degree of resolution possible 
in computing variance and the manner of application to the 
general BKP form. 

Chapters II and III of this study have traced the major 
steps in the development of the three model forms. Although 
some comparative comments have been made, it remains for 
Chapter IV to present the general overview of the models 
and to identify the salient features that might order the 
relative worth of each. Hypothesized examples are proposed 
to allow some objective interpretation of the models that 
has been lacking due to the subjective nature of the material 


presented to this point. 
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IV. COMPARATIVE ANALYSIS AND RESULTS 





The development of both models for determining gun system 
BKP has demonstrated many of the usual obstacles inherent 
with attempts to explain mathematically complex phenomena 
such as military combat operations. The use of assumptions 
to overcome points of resistance or excessive complexity 
serves the function of lubrication in the development, but 
it is always a lubricant with a price. The extent of that 
cost is difficult effect to measure, especially when many 
assumptions are layered within the development. These two 
BKP models are fraught with assumptions that facilitate the 
mathematics of each, but which are at best difficult to 
evaluate with respect to the net effects on the accuracy of 
model results. -It is for this reason that this study makes> 
frequent reference to the need for data from actual system 
tests to assist in validating the models and to measure the 
effects of specific assumptions. Data from tests of existing 
systems may be applied generally to these models without the 
models losing their general applicability which earlier was 
proposed as one of their basic raisons d’etre. This study 
has been conducted without the benefit of test data so that 
assumptions have been necessarily discussed from a theoretical 


and intuitive point of view with conclusions withheld in many 


cases pending validation. 
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Both BKP models are initiated with the realization that 
BKP is one measure of effectiveness that has a broad base of 
application in the field of air defense gun system evaluation 
and simulation. They present and describe identical gun/ 
target environments and engagement procedures whereby one 
gun system encounters one non-maneuvering target. The concepts 
involved in the process of developing both models are essent- 
ially the same as for any general target destruction model. 
The process consists of factors which affect the distribution 
of projectile impact points and the lethality of the projectile 
as a function of impact points and target characteristics. 
The factors which affect the impact point distribution may 
be considered as potential error producers which may be synth- 
esized from sub-components, and the system lethality may be 
expressed as a function of target and projectile parameters. 
These concepts allow the general approach of analysis to 
be one of modeling a random process. The random variables any 
the process are target location error, aim point error, and 
gun system projectile delivery error. In both models, target 
location errors and aim point errors have been combined for 
convenience ~~ error term is assumed to have its distri- 
bution described by some probability density function. Specifi- 
cally, the circular normal and the bivariate normal distri- 
butions have been assumed appropriate for Model's I and II , 
respectively. When firing bursts, the aim point error is 
realized only once for each burst whereas the delivery error 


is realized for every round. This delivery error suggests 
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some correlation between rounds in a burst and some correlation 
between effects of each round when transforming SSKP to BKP. 

Conclusions about a preferred model for BKP would be 
somewhat dangerous and suspect at this point without test 
data to support them, but it is possible to represent trends 
in the general application of the three BKP forms. Sample 
computations have been generated by permuting hypothesized 
values for burst size with varied levels of dispersion for 
the projectile and burst center distributions. Bursts range 
in size from 1 to 100 rounds while vulnerable area is held 
constant at 1.0 square meter. Dispersion, as represented by 
the standard deviations of the distributions discussed, are 
considered in three categories to reflect near, intermediate 
and long range effects. The intermediate range values for 
standard deviation are hypothesized and are then halved and 
doubled to reflect near and long range conditions, respective- 
ly. The selection of actual magnitudes for these values is 
completely arbitrary and in no way reflects real data. The 
range of values selected is intended only to provide a large 
enough spread to be able to identify representative trends in 
the model quantitative results. Tables I, II and III are 
presented to demonstrate the effects on BKP of these permuta- 
tions, and to show trends in the models individually and 
comparatively. 

Results for each of the three BKP forms are remarkably 
similar for burst sizes less than 20 in the intermediate 


range. The rapid increase in values for Model IIA as burst 
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sizes increase past 20 reflect the point at which the assump- 
tion that n x SSKP is small tends to lose its veracity. The 
near and long range cases support the conclusion that "small' 
refers to values of approximately 0.2 or less for n x SSKP. 

The development of Model IIB demonstrated what appeared 
to be some clear improvements theoretically over some of the 
more undesirable assumptions of Model I, but the sample 
results in each of the three tables show near equivalence in 
BKP values throughout. It was suggested earlier that if a 
positive correlation existed between rounds in a burst, then 
the independence assumption in Model I would result in con- 
sistently high estimates for BKP. The result of slightly lower 
BKP values for Model IIB may reflect the degree of affect 
caused by that different approach. 

The point has been emphasized that the models require 
validation before their individual merits can be ascertained, 
however, they do have considerable intuitive and practical 
appeal as they currently exist. The relationship between 
model parameters and solution behavior are easily recognized 
in each model and they are easily applied in simulation. 

They represent a convenient practical and theoretical base 


from which to study or perform parametric analysis on exist- 


ing systems or on systems still on the drawing board. 
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APPENDIX 


Variance of the Critical Random Variables 


Each component of the error vector at time, +t, was assumed 
to be the aggregation of errors attributable to independent 
sources, each of which was said to be normally distributed. 
Consider the development of these random variables. 

i. imrrers Contributing to As. 

a. Azimuth Sensing - Let Al(t) be the random variable 
denoting the angular error resulting from sensor 
azimuth errors. It results in error on the X-axis 


equal to (Al-R)(t). Where, 


E {(Al-R)(t)} = RM,,(t) 
* (GeR Ge) ) = ee on (+) 


b. Azimuth Rate Sensing - Let A(t) be the random 
variable denoting the angular error resulting 
from sensor azimuth rate errors. It results in 
error on the X-axis equal to (A+ to+R)(t) where 
te is the time of flight of the projectile to 
predicted intercept. Mean and variance terms 


are given by: 


E {(A+t,'R)(+)} 


| 
© 


V {(Att,+R)(+)} 


tl 
ct 
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Gun Pointing Error - Let A2(t) be the random 
variable denoting azimuth error resulting from 
gun pointing errors. It results in error on the 
X-axis equal to (A2:R)(t). Mean and variance 


are given by: 
E ((22- Riga) = My o( t) 


v { (A2-R)(t) }= R2. ok o( t) 


Errors Contribumiimne to A EB. 


Errors with respect to the Y-axis which constitute 


AE(t) are developed identically as for , A(t). 


errors Contributing to AR. 


Ae 


Range Sensing - Let Rl be the random variable 
denoting the range error resulting from range 
sensor errors. It results in error on the 2-axix 


equal to (R1)(t). Mean and variance are given 


by: 


E {(R1)(t)} =o (4) 


R1 


V {(R1)(t)} = og (t) 


Range Rate Sensing - Let R be the random vari- 
able denoting range error resulting from range 
rate sensing errors. It results in an error on 
the Z-axis equal to (R+t,) (+t). Mean and variance 


terms are: 
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E ((R+t,)(t) } = tp Mg(t) 
V ((Retp)(t) 3 = +2. 0 F(4) 


c. Muzzle Velocity Error - Let 4 V be the random vari- 
able denoting muzzle velocity variations result- 
ing in range errors ( AR). Reference 1, by 
employing the 3/2 power law, shows that the range 


error contributions may be represented as: 


E {AR } due to AV 


Hi 
© 


V {AR} due to aV = (—R)- of 


Considering these error contributions random variables, it 


may be synthesized that: 





M, = (R'Mg, + RM) (t) 
Mm = (R‘Mp, + R-Mpo)(t) 
i (Mp, + Mg tp) (+) 

And 
a © = (Re. rs; ae ais a te RY. 5 4) (4) 
o = (R295, +R? 102, + 42 «RP Log) (+) 
o® = 55+ ts wine ( v, er 


Once again, with target vulnerability represented ina 


plane, these results are applicable only when they can be 
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Projected into that plane. Combined with Eq. (14) ana Be. 


(15), the distribution parameters become: 


M+ ={R- : gene » © 
eRe OR+ Mes ve ( Mpy + Mp ‘tp, ) }(+) 








‘ = 5 ° _EsR- e e 
M,, {R* Mp, + Re Mpa t ve (Mp, + Mp ‘tp) }(t) 
2 2 2 2 
2 2 
ee * = ‘ : 
{F eee oe +o, oP ) 
2 
A-R 2 e R é 
o ~~ ee ee (—-) Gy)} (+t) 
Pp 
Pp 
2 
D 2 2 2 . 
? ee ee ee 
p 2 2 
m’°R R ) 2 


and the correlation coefficient is given by: 


qi Cov ("A, Aare. 


where 


— 


Corga' AB‘) =E ( (AA -ECAAWUAE - EC AB"'))} 
This is easily expanded to show that: 


= Rp .* _2 
OOGGE A, AE ) = (A°B)( =a) oe 
p 
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such that 


‘es ae 
7 Cee ) 
ue 


ho 


ie 
t 
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These are the input variables which lead to the solution of 
Eq. (17) and the value of SSKP. They demonstrate a technique 
for synthesizing SSKP as a function of its many basic inputs. 
The task of accurately determining and validating each of 
these components remains to be accomplished. 

It 1s encouraging to see such a logical and basic approach 
to the calculation of delivery error variance, but that appeal 
Should not be allowed to cloud the practical aspects of its 
inclusion or adaptation to the model. The actual variance 
values used in the quantification of BKP are engineering 
estimates that may or may not be correct. The question con- 
cerning whether it is better to use one estimated variance 
value for the overall system or a synthesized value accumulated 
from a number of lesser estimates is no less difficult to 
answer in this case then for any other modeling problem. 

The development traced in this appendix appears to have a 
great deal of merit, but it will be of no real consequence 
to the computation of BKP unless it is possible to show 

confidence in the methodology for fixing the values of the 


delivery error variance components. 
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